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ABSTRACT 

This paper studies the long term evolution of planetary systems containing short- 
period planets, including the effects of tidal circularization, secular excitation of eccen- 
tricity by companion planets, and stellar damping. For planetary systems subject to all 
of these effects, analytic solutions (or approximations) are presented for the time evolu- 
tion of the semi-major axes and eccentricities. Secular interactions enhance the inward 
migration and accretion of hot Jupiters, while general relativity tends to act in opposi- 
tion by reducing the effectiveness of the secular perturbations. The analytic solutions 
presented herein allow us to understand these effects over a wide range of parameter 
space and to isolate the effects of general relativity in these planetary systems. 

Subject headings: Stars: Planetary systems 

1. Introduction 

Starting with the discovery of extrasolar planets (Mayor & Queloz f995; Marcy & Butler 1996), 
a substantial fraction of the planetary orbits have been found close to their stars, with periods P ~ 4 
days. These objects are often referred to as "hot Jupiters". With ~ 200 planets detected to date 1 , 
the distribution of orbital periods shows a measurable pile-up at periods P = 3-5 days, i.e., roughly 
10 percent of the currently detected planets have periods P < 5 days. Although the close planets 
are the most easily detected, this finding is not a selection effect: The current statistics indicate 
that 1.2 percent of all FGK stars have hot Jupiters within 0.1 AU of their stars (Marcy et al. 2005). 
These hot Jupiters are subject to tidal interactions with their central stars (e.g., Goldreich & Soter 
1963), and this star-planet coupling can influence the long term evolution of these systems. If the 
system contains additional bodies, then planet-planet interactions can also affect the long term 
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evolution of the close planets. In our companion paper (Adams & Laughlin 2006a; hereafter Paper 
I), we present a treatment of secular interactions for (non-resonant) systems with multiple planets. 
In this paper, we combine this formulation of secular interactions with star-planet interactions to 
study the long term fate and evolution of close planets. 

The basic theory of secular interactions, as applied to extrasolar planetary systems and used 
herein, is reviewed in Paper I (see also Murray & Dermott 1999; hereafter MD99). These interac- 
tions allow for planetary orbits to exchange angular momentum so that orbital eccentricities change 
on secular time scales that are long compared to both the orbit time and observational monitoring 
times (tens of years) , but short compared to the ages of the systems (a few Gyr) . The characteristic 
secular interaction time scales provide a simple metric of the importance of planet-planet pertur- 
bations within a given system. For a collection of observed multi-planet systems, these secular 
time scales fall in the range 100 - 50,000 yr, as listed in Table 1 of Paper I. On longer time scales, 
tidal interactions with the star act to circularize close orbits and can lead to continued inward 
migration. This process might have a bearing on the observed pile-up of planets with periods of 
3-4 days and can lead to considerable energy input into the planetary atmospheres. Although this 
process has been considered previously (e.g., Trilling 2000; Bodenheimer et al. 2001, 2003; Yu &i 
Goldreich 2002; Mardling & Lin 2002, 2004; Gu et al. 2003; Faber et al. 2005), this paper takes 
the additional step of providing analytic expressions for the evolution time for various classes of 
systems, including both tidal damping and secular interactions. This work thus provides additional 
analytic insight into the problem as well as elucidating the role played by secular interactions. The 
version of secular interaction theory formulated in Paper I includes the leading order corrections 
for general relativity (GR), which causes the periasta of planets to precess forward in their orbits. 
As a result, this semi- analytical treatment also allows us to explicitly delineate the role played by 
general relativity in the long term evolution of close planets. 

For a single planet in a close orbit, the tidal circularization effect is usually written in terms 
of the time scale for eccentricity decay, r c j r = —e/e. For close orbits in multiple planet systems, 
however, the eccentricity is excited via interactions that we model according to secular theory. The 
first approximation to the long term behavior is to assume that the orbit decays with constant 
angular momentum and that the semi-major axis decreases on the same time scale, i.e., 

a_ 2ee 2e 2 1 

a ~ ~ ~ (T^2) — > u 

where the time scale T c i r for circularization can be written in the form 

where Qp ~ 10 5 — 10 6 is the tidal quality factor and Rp is the planet radius. For extrasolar planets, 
the quality factor Qp and the radius Rp depend sensitively on the planetary mass, temperature, 
and composition (Bodenheimer et al. 2003). This general form for the time circularization scale is 
well known (e.g., Goldreich & Soter 1966), but includes additional factors to account for the effects 
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of nonzero eccentricity (see Hut 1981). The choice F « 1 + 6e 2 + C(e 4 ) provides a reasonable 
approximation to the results for close binaries (Hut 1981) for moderate eccentricities. Inserting 
representative values, and taking the limit e — > 0, we write the circularization time scale in the 
form 

(&) a ©~ 3/2 (IT (^r ■ « 3 » 

Throughout this paper we work in terms of the dimensionless time scale 

t = t/r . (4) 

These systems are also affected by a stellar damping effect in which energy is dissipated in 
the star due to tides raised by the planet. The effectiveness of this process is determined by an 
analogous parameter Q*, the tidal quality factor of the star. The net result of this process is to 
cause the semi-major axis of the inner planet to decay on a time scale Tp* (again, see Goldreich & 
Soter 1966, Hut 1981) that can be written in the form 

— ^ r. ?(*)•(£) (==)'. 

As defined here, T is generally a small parameter. For Jupiter-like planets, Rp ~ 0.1-R* and 
mp ~ 10 _3 M* so that Tj ~ 0.03(Qp/Q*). For Neptune-like planets, we find Tj^ ep ~ 0.015(Qp/Q*), 
but we expect the Qp value to be much smaller so that T^ ep <C Tj. 

This paper considers hot Jupiters subject to tidal circularization effects (with strength deter- 
mined by Qp), both in single planet systems (§2.1) and in two-planet systems where the eccentricity 
of the inner planet is excited through secular interactions (§2.2). Next we consider hot Jupiter sys- 
tems with additional stellar damping effects (with strength determined by Q*), both for single 
planet systems (§3.1) and two-planet systems (§3.2). Our results are summarized in §4 along with 
a discussion of their ramifications. 



2. Long Term Evolution with No Stellar Damping 

The basic goal of this section is to provide an analytic understanding of the long term behavior 
of the hot Jupiter systems in the absence of stellar damping terms (note that stellar damping is 
included in the following section). For single planet systems that experience tidal forces only, we 
define g(t) = ai(i)/ai(0). Similarly, for two planet systems, we define f(t) = ai(t)/a\(0), where 
the inner planet experiences circularization from the central star and eccentricity excitation from 
the other planet. 
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2.1. One planet system with no stellar damping 

For a one planet system, the evolution is described by two coupled equations of motion. If we 
work to the same order of approximation in e 2 , these equations take the form 

dg „ 2 -11/2 , de _ 13 / 2 , . 

— = -2e 5 ' and -j = -eg 7 , (6) 
at at 

which be be combined to form the second order differential equation 

dg l( fg VLdg = 13/2 

W dt 2 2gdt g ' { } 

The first integral of this equation can also be found, i.e., 

| = -2 5 - 11 /2 (e 2 (o)+ln5)) (8) 

where ei( ) is the eccentricity of the inner planet at t = 0. In contrast to the two planet case, 
where continued eccentricity forcing by the second planet leads to continued evolution, this system 
asymptotically approaches a minimum value of semi-major axis, i.e., a minimum value of g given 
by ^oo = ex p[ — e i(o)]- Combining the two expressions in equation (6), we can also solve directly for 
the eccentricity as a function of the factor g, i.e., e 2 = e 2 ^ + In g. This expression does not satisfy 
conservation of angular momentum exactly because of the approximation made at the start (where 
we work to only leading order in e 2 ). 

If we include the additional factors to enforce conservation of angular momentum, the equations 
of motion take the form 

and the resulting second order differential equation becomes 

(%- 1 ^ + 1 ^ = -2g-^^- 2 . (10) 
V dt' dP 2g dt y 1-4 K ' 

Here we have eliminated the eccentricity dependence using conservation of angular momentum, i.e., 
g(l — e 2 ) = (1 — Cq) = J 2 = constant. This second order differential equation can be integrated 
once to find the time as a function of g, i.e., 

2 J g g-.P ■ 

The integral can also be evaluated to find the implicit solution 

t = £(1 - S 11 ' 2 ) + £(1 - 9 9/2 ) + y (1 - 9 712 ) + £(1 - / /2 ) + 

£(1 - ^) + ^(1 - 9 '/>, + ij» in • ^±£| . (12) 

Both the time integral (eq. [11]) and the solution (eq. [12]) indicate that the planet only reaches 
its final location (given by g = J 2 ) asymptotically in time (t — > oo). 
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2.2. Two planet system with no stellar damping 

For the long term evolution of a two planet system, the eccentricity of the inner planet is 
driven by the interactions with the outer planet. Within this set of approximations, the evolution 
is given by the time-averaged equation 

^ = -2(e 2 )(— ) = -2<e 2 > r 13/2 , (13) 

where the angular brackets denote time averages over an intermediate time scale that is long 
compared to the secular time scales (e.g., see Table 1 of Paper I) and short compared to the 
circularization time scales (eq. [3]). As the system evolves, the semi-major axis of the inner planet 
decreases, which in turn causes the secular averaged square eccentricity (e 2 ) to decrease. Notice 
that by time-averaging the square of the eccentricity and the circularization time scale as separate 
quantities, we are making an approximation that limits the accuracy of this treatment to 0(e 2 ). 
In practice, however, we evaluate r c j r in the e — ► limit, and the secular interaction theory used 
here is only accurate to second order in e (Paper I; MD99), so that this order of approximation is 
consistent with our general framework. 

In the absence of other effects, the equation of motion for the two planet system can be written 
in the form 

| = - Vr 11/2 , (14) 
at 

where the eccentricity excitation amplitude r\ takes the approximate form 

16 (1 + n r 3 - 5V^f) 2 + (25/4)<54 /2 /5/2 ' 

where we have defined 5 = m\jrwi and ao = ai(0)/a2(0). The dimensionless parameter IIo = 
AGMla\j (rri2C 2 a\) provides a measure of the importance of general relativity in secular interactions, 
where the semi-major axes are evaluated at t = (from eq. [17] of Paper I; see also Adams & 
Laughlin 2006b). The quantity e 2 (o) is the eccentricity of the second planet evaluated at t = 
(although this eccentricity is assumed constant in this set of approximations) . It is useful to collect 
the constants into the dimensionless composites 

^ = ~^ e 2(o) a o > B = 5^/a^, and C = ^-5a^ 2 . (16) 



With these definitions, the equation of motion can be written in the form 

df Afl* 



(17) 



di (no + P - Bpiiy + c/ 17 /2 • 

By absorbing the leading coefficient A into the time variable, we can define the relevant time scale 
ts for the evolution of this system, i.e., 

T5 = 25a 2 e 2 = A ' (18) 
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This time scale is often much longer than the circularization time tq that applies for single planet 
systems. For example, using the observed parameters for the inner two planets of the Ups And 
system, we find A « 0.0022 so that ts/tq rs 454. 

The equation of motion (17) can be integrated to provide the solution in implicit form 
At = 1(1- f' 2 ) + |Do(l - / 3/2 ) + \^{r" 2 ~ 1) + j (1 - f) 

+^B\l - f 11 ' 2 ) - U B(1 - f 2 ) - 2 -B(l - f) . (19) 
In the limit B, C <C 1 (which often holds), this solution can be simplified further to the form: 

At = 2 -(l- f' 2 ) + |Do(l - f' 2 ) + 2 -Ul(f-^ - 1) . (20) 

From this expression, we can read off the asymptotic behavior. As the time t increases, /<1, and 
the third term in the equation dominates. As a result, as t — > oo, the function f(t) approaches the 
limiting form 

f(t) -+ < 3 (3Ai/2)- 2 / 3 . (21) 

Unlike the case of a single planet system (see §2.1), the migrating planet can continue to lose energy 
and can become arbitrarily close to the origin. In practice, however, once the planet reaches the 
stellar surface where / = / m j„ = i?*/ai(0), the planet will be destroyed and evolution will be over. 
The time required for the planet to reach the stellar surface is thus given by 

AU = 2 -(l - fJl) + |no(l - fJl) + ^(/-;f - 1) . (22) 

Notice the important role played by general relativity in this setting. In the absence of relativistic 
corrections, the (leading order) solution would have the form f(t) = (1 — 9At/2) 2 ^ 9 , which reaches 
/ = in the relatively short time t* = 2/9 A. Relativistic precession thus acts to keep the planet 
from being accreted by the star. This claim can be quantified by inserting typical values of interest; 
let ai(0) = 0.05 AU and R* = 1.0 Rq so that the total evolution time given by equation (22) can 
be written 9At*/2 «1 + 6Ho + 105U"q. With these values, e.g., the condition for relativistic effects 
to dominate becomes n > 0.073. 

The portion of parameter space for which two planet systems lead to significant (short) accre- 
tion times is depicted in Figure 1. In this application, we assume that the inner planet is a hot 
Jupiter, with mass mi = mj, radius Rp = Rj, and starting semi-major axis a\ = 0.05 AU. The 
eccentricity of the inner planet cycles through e\ = 0, although (e) 7^ due to secular interactions 
with the second planet. The stellar mass M* = 1.0 Mq and the tidal quality factor Qp = 10 6 so 
that the circularization time scale tq = 1.6 Gyr (see eq. [3]). The figure shows the orbital elements 
of the second planet (02,62) required to drive the hot Jupiter into the stellar photosphere over a 
fiducial time scale of 5 Gyr. The accretion time depends on the mass of the second planet, which 
is assumed to vary over the range 777-2 = 1 — 5?77j, corresponding to the five solid curves shown in 
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Figure 1. The region of the a — e plane above the curves represents solar systems in which the inner 
planet would be driven into the star on time scales less than 5 Gyr. The dashed curve delineates 
the (much larger) region of parameter space for which the inner planet would be driven into the 
star in the absence of general relativity (where mi = raj). Comparison of the dashed curve with 
the upper solid curve thus illustrates how GR acts to prevent the accretion of planets. 



3. Long Term Evolution with Stellar Damping Term 

Given the solutions discussed above, the next correction is to consider the dissipation of energy 
that occurs due to the planet raising tides on the star. This section presents solutions for the long 
term evolution of these systems, including both tidal circular ization (due to energy dissipated in 
the planet via Qp) and orbital damping (due to energy dissipated in the star via Q*)- We present 
solutions for single planets systems (§3.1) and for two planet systems in which the bodies are also 
subject to secular interactions (§3.2). 



3.1. One planet system with stellar damping term 



In the presence of stellar damping, the equations of motion for the one planet system take the 



form 



dg 

-4 = -2— 
dt 1 



7 -n/2 _ rv 11 /? 



and * = -eg- 13 / 2 
dt 



(23) 



In the limit of zero eccentricity, the stellar damping term acts to decrease the semi-major axis of 
the planet according to g(t) = (1 — 13rt/2) 2 / 13 ; the system thus has a "natural" damping time 
scale of t c = 2/13r. Combining the above equations, we can solve directly for g(e), 

Hirfc!)- <-) 

We can insert this form back into the eccentricity evolution equation and solve for the time as a 
function of e, 



T 



2 \ 13/2 



In the limit of small eccentricity, we can evaluate the integral to find the approximate solution 

2 



u (1 - 4)13/2 w I 1 ~ (e/eo)13r/2 

Alternatively, we can find the approximate solution for e(t), 



e(t)«eo[l-yr(l-eg) 



-13/2, 



2/i3r 



(26) 



(27) 



This function, in conjunction with the solution for g(e) found above, thus specifies the semi-major 
axis of the planet as a function of time. 
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The total evolution time is given by the integral of equation (25) in the limit that e — ► 0. This 
time scale can be written in terms of the series 

(T . (1 _ e ^ E ^_^_, (28) 

n=0 ' 

where the coefficients b n are defined through the recursion relation 

b n+ i = (13r/2 + n) b n with b = 1 . (29) 

These analytic results tell us a lot about the long term behavior of the system. Equation (24) 
implies that the planet does not reach the origin (g — > 0) until its eccentricity vanishes. Further, 
the integral in equation (25) is finite for all T > 0, so the evolution takes place over a finite time 
(equivalently, the series in equation [28] converges). Finally, equation (28) provides a good working 
approximation to the total evolutionary time scale, i.e., 



13 



UTel 13r(2 + 13r)eg 



+ -: tt^ + 



i3r 4 + i3r 4(8 + i3r) 



(30) 



For many applications, the first term provides an adequate approximation. Notice that the total 
evolution time is shorter than the characteristic time scale rp* = tq/T (i.e., t = 1/T), as defined by 
equation (5). When the leading order term is a valid approximation, the evolution time is shorter 
than the characteristic time t c = 2/13r by an additional factor of J 13 = (1 — e^) 13 / 2 . For example, 
if eo = 0.28 (the median of the observed sample of extrasolar planets), the eccentricity causes the 
evolutionary time to be shorter by a factor of ~ 2 compared to evolution with no initial eccentricity. 

The observed sample of extrasolar planets shows a population of planets with periods of 3 - 4 
days, but a deficit of planets with shorter periods. One could, in principle, explain this signature 
if the shorter period planets were all accreted by their central stars. If this explanation were true, 
then the total evolution time defined above must be comparable to the stellar ages r* (which are 
typically several Gyr). The requirement that all planets with (initial) semi-major axis less than 
ao are accreted within the stellar age r* can be evaluated by using the leading order term for the 
evolution time (in eq. [30]) and the definitions of To and T (eqs. [2] and [5]). The result takes the 
form 

where we have assumed solar properties for the star (M* = 1M and i?* = 1R&). Thus, planets of 
roughly Jovian mass can be accreted in a typical stellar age (5 Gyr) provided that ~ 10 6 (which 
is a reasonable value, implied by considerations of eccentricity damping in close binaries, e.g., Hut 
1981). Because of the sensitive dependence on semi-major axis a, planets that are accreted spend 
relatively little time with shorter periods, and, the predicted cutoff in observed period is relatively 

13/2 13/3 

sharp. Note that r* oc a oc P Q 1 , so that a planet with a 2 day period is accreted 20 times faster 
than a planet with a 4 day period. For a given age of the stellar population, one would expect to 
see far more planets in 4 day orbits than in 2 day orbits (if no observational biases are present). 
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When stellar damping is included, some of the orbital energy is dissipated in the planet, and 
some is dissipated in the star. For purposes of finding the impact of energy dissipation on the 
planet (for example, the effect on the planetary radius), we need to find the fraction of the energy 
that is dissipated in the planet. For this case, we can consider both the energy dissipation and the 
evolution of semi-major axis to be functions of the eccentricity. Equation (24) specifies the function 
g(e). The energy dissipated in the planet, as a function of eccentricity e, can be written in the form 



AE P GM*mp 



oo(l-e§)(2-r) 



e 



(e /e) r l . (32) 



This function of eccentricity can be converted into a function of time using the time evolution 
equations (25 - 31) discussed above. The self-gravitational energy of a Jovian planet can be 
written Eqp = i]Gmp/ Rp, where the dimensionless parameter rj ~ 5/3 (although the exact value 
depends on the internal structure and of the planet). We can thus determine the conditions for 
which a planet will experience a dissipational energy AEp that is comparable to its self-gravity 
Egp- 



AE P M* R P 



We) 



E GP m P a V (l - eg) (2 - T) 

In the second approximate equality, we have assumed M* = 1.0M & and Jovian properties for the 
planet (mp = 1.0m j and Rp = l.ORj). If a hot Jupiter starts with ao = 0.05 AU and eccentricity 
eo ~ 0.51, then the amount of energy dissipated within the planet through tidal interactions is 
comparable to its self-gravity, so that substantial structural changes can be forced upon the planet. 
In most cases, however, we expect a smaller eccentricity, so that only a fraction of the planet's 
self-gravitational energy would be dissipated (as given by eq. [33]). For the median of observed 
eccentricity of the current sample, eu ~ 0.28, the ratio AEp /Egp ~ 0.24, which is still large 
enough to be significant. 



3.2. Two planet system with stellar damping term 

For a two planet system with secular interactions and stellar damping effects, the equation of 
motion takes the form 

df _ Af 5/2 w -n/2 



di (no + / 3 - Bf 7 / 2 ) 2 + cf 17 / 2 



Yf- lv ' 2 . (34) 



If we make the same approximations as before (§2.2), where B, C <C 1, the equation of motion can 
be written in terms of a formal implicit solution of the form 

Although the integral cannot be evaluated in terms of elementary functions, a good working ap- 
proximation can be found if we assume that the first, secular term in the equation of motion is 
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important only for / ~ 1. To find this approximation, we first rewrite the integral in the form 



i; 



Tt= / df 



f n/2 _ (A/T)! 27 ' 2 

1 (n + /3)2 + ( A/r)/8 



(36) 



We then assume that the second term (which arises from secular perturbations) only is important 
at the beginning of the evolution with / ~ 1, so we can set / = 1 in the denominator. The resulting 
expression for the evolution time becomes 

This approximation is valid as long as the ratio A/T is not too large. For example, the error is less 
than 15% as long as A/T < 4 for Ho = 1. In the extreme limit r < 1, where this approximation 
fails, the additional damping term does not play a role and we can use the no damping solution 
as a good working approximation. Comparison of this result with that of the previous subsection 
indicates that in order for the secular interaction of the second planet to play a role, the second 
term must compete with the first. Almost equivalently, the first term in the evolution equation 
(34) must compete with the second, i.e., A ~ T(l + ITo) 2 . Inserting typical masses and radii for 
the planets and star, the requirement for secular interactions to dominate the long term evolution 
becomes 

Since a\ <C 02 for typical systems, this requirement is not met unless Q* » Qp. As a result, the 
stellar damping effects tend to dominate secular interactions as the inner planet migrates closer to 
the central star. 

For solar systems with sufficiently small semi-major axes and high eccentricities, secular in- 
teractions coupled with stellar damping drive the inner planet into the star. To define a region of 
parameter space for which this effect is important, we take the inner planet to be a hot Jupiter 
with mi = mj, a = 0.05 AU, and an eccentricity that cycles through e\ = 0. The orbital elements 
for the second planet that lead to planetary accretion within 5 Gyr is shown in Figure 2, where we 
have taken = 10 6 . Since the effect depends on the mass of the second planet, we show curves for 
m2 = 1, 3, and 5 mj. This plot is thus the analog of that shown in Figure 1, which does not include 
the stellar damping term. As expected, inclusion of stellar damping effects leads to a larger region 
of parameter space for which planets can be accreted. The dashed curve delineates the portion of 
parameter space for which a system with rri2 = mj would drive its inner planet into the star within 
5 Gyr in the absence of relativity. 



4. Conclusion 



This paper presents solutions for the the long term evolution of four types of planetary systems 
- one and two planet systems with tidal circularization, both with and without the inclusion of 
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stellar damping effects. For the two planet systems, planet-planet interactions are modelled using 
secular interaction theory (MD99) including leading order relativistic corrections (Paper I). The 
solutions are presented in implicit form; for example, t(f) gives the time as a function of the inner 
planet's semi-major axis, as specified by / = ai(t)/ai(0). In the absence of stellar damping, we 
find exact analytic solutions to the long term evolution, as given by equations (12) and (19). For 
systems in which the stellar damping term is important, we find approximate analytic solutions 
given by equations (26) and (37). 

In multiple planet systems, secular interactions enhance the accretion of inner planets by the 
central star, provided that the outer planet has a sufficiently large mass, small semi-major axis, 
and/or large eccentricity. We have presented a quantitative assessment of the orbital elements of 
a second planet required to drive a hot Jupiter into the central star within 5 Gyr, for two-planet 
systems both with (Fig. 2) and without (Fig. 1) stellar damping. General relativity acts to delay 
the accretion of planets, i.e., the region of parameter space for which the inner planet would be 
accreted within a given time would be much larger in the absence of relativistic effects (see the 
dashed curves in Figs. 1 and 2). 

For multiple planet systems with hot Jupiters as inner planets, secular interactions tend to 
give the inner planet nonzero eccentricity. This continual addition of eccentricity, in conjunction 
with the circularization processes experienced by such planets, leads to large amounts of energy 
dissipation within the hot Jupiters (eqs. [32] and [33]). Such extreme dissipation, in turn, may 
lead to mass loss and thereby explain the relatively small masses observed for close planetary 
companions. This process should be studied in greater detail in the future. 

The framework of analysis presented here can be readily applied to individual systems whenever 
the orbital elements of the constituent planets are known to reasonable accuracy. We maintain an 
up-to-date catalog of the known extrasolar planets 2 where the time scales derived here are tabulated 
using the best available fits to the radial velocity data sets. This data base, in conjunction with 
the analytic solutions presented here, should provide a useful resource for further research on close 
planetary systems. 
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separately from Paper I. This work was supported at U. Michigan (FCA) by the Michigan Center for 
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and was also supported at U. C. Santa Cruz (GL) by NASA through the Terrestrial Planet Finder 
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Fig. 1. — For systems that contain a hot Jupiter, this plot shows the orbital elements of the second 
planet (02, e-i) required to drive the hot Jupiter into the stellar photosphere over a time scale of 5 
Gyr. All cases assume that the inner planet has mass m\ = mj, semi-major axis a\ = 0.05 AU, 
and an eccentricity that cycles through e\ = (where (e) / due to secular interactions with the 
second planet). The five solid curves correspond to five choices of the mass of the second planet, 
as labeled. The region of parameter space above the curves corresponds to systems in which the 
inner planet is driven into the star on shorter time scales t < 5 Gyr. The dashed curve delineates 
the (much larger) region of parameter space for which the inner planet would be driven into the 
star in the absence of general relativity (for 771,2 = fnj). 
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Fig. 2. — For systems that contain a hot Jupiter, this plot shows the orbital elements of the second 
planet (02, e^) required to drive the inner planet into the star within 5 Gyr, where a stellar damping 
term is included with Q* = 10 6 (compare with Fig. 1). The inner planet has mass mi = mj, semi- 
major axis ai = 0.05 AU, and an eccentricity that cycles through e\ = (where (e) / due to 
secular interactions with the second planet). The three solid curves correspond to different masses 
of the second planet, as labeled. The region of parameter space above the curves corresponds to 
systems in which the inner planet is driven into the star on time scales t < 5 Gyr. The dashed 
curve delineates the region of parameter space for which the inner planet would be driven into the 
star in the absence of general relativity (for 777-2 = mj). 



